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Abstract 

We discuss energy dependence of rapidity gap survival probability in the 
double-pomeron exchange processes with account of the reflective scatter- 
ing effects. 
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Introduction 



Double-pomeron exchange processes are the subject of the active theoretical stud- 
ies since the beginning of the nineties, when the possibility to discover the Higgs 
production in these processes was discussed in the seminal papers 03I21I3. Nowa- 
days these studies are motivated by the current experiments at Tevatron [4] and 
planning for the future experimental studies at the LHC [5J. Needless to say that 
the processes with large rapidity gaps are important tool in the search of new 
physics. The extensive studies of the magnitude and energy dependence of gap 
survival probability have been performed, results of these studies together with 
the estimates of the Higgs production cross-section can be found, e.g. in papers 
IIS H El EH QH. Recent results and review can be found in E2. 

Perturbative QCD calculations can be used for calculations of the cross-sections 
of such processes where hard scale such as a large mass is present. But the ne- 
cessity to take into account soft rescatterings in initial and final states leads to 
uncertainties related to the absence of the theoretically well established proce- 
dure of the multiparticle amplitude unitarization which is commonly reduced to 
the absorptive corrections due to the rescattering processes. Dynamics of such 
interactions in hard diffraction processes is accounted then by introduction of a 
factor which is known as a gap survival probability [0, i.e. a probability that two 
colliding particles pass through without any interaction. 

The gap survival probability (SP) was denoted as (\S\ 2 ) since it was defined 
by the following relation : 

{lb{) ~ / °° D H (b)dPb ' ( ) 

where Dn(b) is the probability to observe a specific hard interaction in collision 
of the hadrons hi and h 2 and \S(s,b)\ 2 , where \S(s,b)\ 2 = exp[— 0(s, &)], is a 
probability that colliding particles do not have any interaction. The above def- 
inition of SP proceeds from the eikonal scattering picture |[T3ll which is usually 
used for the estimates of (|5| 2 ), where the eikonal function is denoted by £l(s, b). 
It implies that the soft scattering picture at very high energy is an absorptive one 
(S > 0) and the elastic scattering amplitude in the impact parameter representa- 
tion does not exceed black disk limit which is a half of the unitarity limit. It is 
worth to note that the diffractive processes (with exclusion of elastic scattering) 
are inelastic processes themselves and therefore give contribution to absorption. 
Account for the soft processes is model dependent and strongly depends on the 
details of the particular model. Estimates for the gap survival probabilities can 
vary in the wide range, e.g. from 0.4 to 22 % [fT4l . Of course, such wide variation 
of theoretical predictions is not encouraging for the experimental studies. 
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In the eikonal picture S(s, b) can easily be treated as the elastic scattering S- 
matrix, however there are theoretical problems with the definition given by Eq. 
(0Q) and with such interpretation of the function S if one goes beyond the eikonal 
scattering picture . As it was shown in |[T5l . at very high energies the elastic scat- 
tering can go beyond the black disk limit |[T6ll and become a reflective at small 
impact parameters, i.e. elastic S matrix can have negative values. The generic 
geometrical picture at fixed energy beyond the black disc limit can be described 
in this case as a scattering off the partially reflective and partially absorptive disk 
surrounded by the black ring which becomes grey at larger values of the impact 
parameter. The evolution with energy is characterized by increasing albedo due 
to the interrelated increase of reflection and decrease of absorption at small im- 
pact parameters. Asymptotically, picture of particle collisions with small impact 
parameters looks like collisions of hard balls. This picture implies that the scatter- 
ing amplitude at the LHC energies is beyond the black disk limit at small impact 
parameters (elastic S'-matrix is negative) and it provides explanation for the reg- 
ularities observed in cosmic rays studies lfT5ll . The straightforward use of Eq. 
ffl l[T7l in this case leads to the non-monotonic energy dependence of SP and its 
increasing dependence. It approaches unity at very high energies which has no 
physical meaning. This limit is merely a manifestation of the pure reflective elas- 
tic scattering (S — > — 1) at s — > oo and 6 = 0. The reflective elastic scattering does 
not result from absorption and therefore the definition given by Eq. (Q~|) cannot be 
applied for the S'-matrices having reflective scattering mode. The limitation of the 
definition Eq. (Q~|) to absorptive scattering is absent if one uses the more general 
definition for the SP (we denote this quantity as R? s to differentiate it from the Eq. 
©) given in (HOSED 

<7di//(hard) 

where <Jdiff (hard) can, in principle, be calculated with perturbative QCD, while 
(Jdiff (phys) is the real cross-section, which can be confronted with the cross- 
section measured experimentally. The definitions Eqs. (Q]) and © are equivalent if 
S is nonnegative in the whole range of variation of the variables s and b. However, 
they are not equivalent if S can become negative at some values of the energy and 
impact parameter. 

We use therefore the definition Eq. © and calculate the energy dependence 
of R%(s) in the unitarization scheme which includes reflective scattering in the 
double-pomeron exchange processes (DPE). 

Now a few comments are to be added. First of all, of course, the quantity R? s is 
model and process dependent. Another important point was mentioned in [9], this 
quantity is probability of surviving rapidity gaps, it is not the probability of pro- 
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duction and surviving rapidity gaps. It is useful to have it in mind under discussion 
of the available experimental data on hard hadronic diffraction at Tevatron. 

Experimental studies at Tevatron energies y/s = 630 GeV and ^fs = 1800 
GeV show that in the hard double-diffractive (DD) interactions the ratio of the 
of the corresponding event numbers decreases with energy [4J. The studies of 
inclusive DD at CDF [HI leads to conclusions on the universality of SP in soft 
and hard diffractive processes and on the absence of suppression of the additional 
gaps compare to the diffractive processes with single rapidity gap, i.e. SP does 
not depend on number of gaps. 

SP in double-pomeron exchange processes; effects of 
the reflective scattering 

As it was already mentioned, almost all estimates of SP are performed on the 
basis of the eikonal amplitude unitarization scheme. To strengthen the obtained 
predictions it is desirable to have predictions based on the other forms of the 
amplitude unitarization. It is important to move this way in the situation when 
there is no unique scheme for the unitarity account. 

Alternative form of unitarization uses the rational representation for the scat- 
tering amplitude. The rational form of unitarization in quantum field theory is 
based on the relativistic generalization EOl of the Heitler equation [21 J. In this 
U -matrix approach the elastic scattering matrix in the impact parameter represen- 
tation has the form: 

l + iU(s,b) 

S{s,b) = - , (3) 

1 — iU{s, b) 

where S(s, b) = l+2z/(s, b) and U (s, b) is the generalized reaction matrix, which 
is considered to be an input dynamical quantity similar to the eikonal function. 
Unitarity equation rewritten at high energies for the elastic amplitude f(s, b) has 
the form 

Im/(s, b) = \f(s,b)\ 2 + rj(s,b) (4) 
where the inelastic overlap function 

is the sum of all inelastic channel contributions. Inelastic overlap function is re- 
lated to U (s, b) according to Eqs. © and © as follows 

Vis b) = lmU{3 > b) (5) 
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It should be noted that 

lmU(s,b)=J2Un(s,b), (7) 

n>3 

where U n (s, b) is a Fourier-Bessel transform of the function 

1 r n ri 3 a - 

Un(s,t) = / JT_^^(Vg i -p a -p 6 )C/:( gi ,....,g n ;p' a ,^)- (8) 
mi Qjn 

J j=l i=l 

U n (qi,....,q n ;p a ,p b ). 

Here the functions C/ n (gi, q n ;p a ,Pb) and q n ;Pa>,Pb>) correspond to 

the ununitarized (input or "Born") amplitudes of the process 

a + b — > 1 + .... + n, 

and the process with the same final state and the initial state with different mo- 
menta p' a and p' b 

a' + b' -> 1 + .... +n, 

respectively. They are the analogs of the elastic [/-matrix for the processes 2 — ► n. 
The sum in the right hand side of the Eq. © runs over all inelastic final states 
|n) which include diffractive as well as non-diffractive ones. Graphically, these 
relations are illustrated in Fig. 1. 



Im U 




n>3 



Figure 1: Unitarity for the function ImU . 
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Then the inclusive cross-section of the process ab — ► cX has the following 
form ll22l . which is similar to the expression for the total inelastic cross-section 
©: 

da _ f 50 . „ I(s,b,q) 



E- 



S7T 



bdb- 



12' 



(9) 



A ""Jo """I 1 -^M)l 

where J(s, 6, g) is the Fourier-Bessel transform of the functions which are defined 
similar to Eq. © but with the fixed momentum q and energy E of the particle 
c in the final state (Fig. 2). It should be noted that the impact parameter b is 




Ks,t,q) 

q\ /q 




n>3 



Figure 2: Unitarity for the function I(s, t, q). 



the variable conjugated to the variable yf—i, where t = (p a — p' a ) 2 and that the 
following sum rule is valid for the I(s, b, q) 

J ^-I(s, b, q) = n{s, b)lmU(s, b). (10) 

The impact parameter b is related to the impact parameters of the secondary par- 
ticles by relation [|23l 

n 

b = j2^ (id 

i=l 

where Xi stands for the Feynman variable x of the i-th particle. We consider the 
double-pomeron exchange processes of the following type 

pp^p + X + p, (12) 
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where + denotes gap in rapidity and X can be a particle or a system of particles. 
To obtain the cross-section of the DPE process (TT2l) we should single out in the 
sum over all inelastic final states \n) the specific final states \n)d pe , which are the 
final states corresponding to the specific kinematics of the corresponding DPE 
process. Then denoting via ( a set of kinematical variables which characterizes 
the kinematics of the final state in reaction (TT2l) . a similar to Eq. © expression 
can be written for the cross-section da dpe /d( 

r^ ,Mt 6 'fl - d3) 



d( J \l-iU{s,b)\' 

The formula (PT3l is consistent with the experimental result of CDF on the inde- 
pendence of SP on the gap number. The function U(s, b) is related to the elastic 
scattering amplitude and can, in principle, be extracted from the elastic scatter- 
ing experimental data, however, in practice, the various models can be used to 
construct the function U(s,b) and get its form at the energies where the data are 
absent at the moment. The models also needed to construct Id pe (s, b, (). The ratio 

J dCI dpe (s, b, 0/ J ^fl(s, b, q) = c(s, b) (14) 



should be considered as a gap production probability. From the Eqs. (flOl) and CH 
it follows that 



/ 



d(I dpe (s,b,() = c(s,b)n(s,b)lmU(s,b). (15) 



To get quantitative results we use earlier constructed model for the U -matrix 
based on the ideas of chiral quark models ||24||. The picture of a hadron consisting 
of constituent quarks in the central part of hadron embedded into quark conden- 
sate implies that overlapping and interaction of peripheral clouds occur at the first 
stage of hadron interaction. The impact parameter dependence of this interaction 
is determined by the function Dc(b), which is a convolution of the two conden- 
sate distributions D^(b) and D b c {b) inside the hadrons a and b. We consider for 
simplicity the case of a pure imaginary amplitude, i.e. U — > iu. The function 
u(s,b) is represented in the model as a product of the averaged quark amplitudes 

N 
i=l 

in accordance with the assumed quasi-independent nature of the valence quark 
scattering, N is the total number of valence quarks in the colliding hadrons. 
The essential point here is the rise with energy of the number of the scatterers 



{s,b) = H(f Qi (s,b)) (16) 
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like \/s. The 6-dependence of the function (fg) has a simple form (/q(6)) oc 
exp(— rriQb/£). The generalized reaction matrix gets the following form 

U (s,b) = g(l + a^-) exp(-^f), (17) 

where M = J2q=i m Q- Here niq is the mass of constituent quark, which is 
taken to be 0.35 The model provides linear dependence on \fs for the total 

cross-sections, i.e. a tot = a + Cy/s in the limited energy range y/s < 0.5 TeV in 
agreement with experimental data Il25ll26ll and asymptotically, o tot oc In 2 s, while 

(Jinel OC Ins. 

Inelastic diffractive processes in the central region (DPE) occur in the periph- 
eral region of the impact parameter space. Constituent quarks, which are respon- 
sible for elastic scattering, are supposed to be located in the central part. This fact 
and relations (TT4)) and (fT5l) allow one to assume factorization at the U -matrix level 
and write down the expression for Id pe (s, b, () in the form 

Id P e(s,b,() = a d p e (s,b,()ImU(s 7 b), (18) 

where ad pe (s, b, () describes hard or soft "Born" cross-section of DPE processes. 
It should be stressed that the factorization ([TBI is assumed to be valid for both 
soft and hard DPE processes at the [/-matrix level. We do not expect similar 
factorization, e.g. for DD processes. This factorization is different also from the 
factorization discussed in [fP9l at the S'-matrix level which is considered to be valid 
for all diffractive processes with hard scale. 

We consider hard DPE processes because they are the most interesting ones, 

i.e. 

(?dpe(s,b,() -> cr H (s,b,()- 

Thus, under construction of cr H (s, b, £) we will suppose that the interaction occur 
at small distances. We can write down then the probability of hard interactions in 
the model as a convolution 

a H (s,b,() = J D^ibJaUs^ + bx-b^OD^ib^dbz, (19) 

where cr%(s, b + bi — b 2 , C) is the cross-section of the hard condensate (parton) 
interactions. If the interaction is hard, than it is natural to choose 

a° H (s, b + bx - b 2 , C) ^ a° H (s, C)S 2 (b + b x - b 2 ). 
'Other parameters have the following values: g — 0.24, £ = 2.5, a — 0.56 • 10 -4 . 
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It means that hard interaction cross-section can be represented as a convolution 
integral 

a H (s,b,()^(T H (s,0 f ^(b-bO^CbOdbi, (20) 

which is similar to the model originally proposed by Chou and Yang lfT3ll . This 
model can be considered, in fact, as a model based on hard parton interactions. 
It is useful to note here that the hardness of the inelastic interaction does not 
imply the smallness of the impact parameter b (cf. Eq. CCD)), the typical values 
of the impact parameter b are determined by the condensate distributions D^ 1 
and D^ 2 which are responsible for the central diffractive production and have 
mostly peripheral impact parameter dependence. The two-scale picture of hadron 
scattering proposed in |fT9ll27l . where the core of hadron is responsible for hard 
interactions and peripheral part - for soft interaction, should definitely take place 
in elastic diffraction, while the hard interaction in the inelastic diffraction can 
take place in the peripheral region as well. Of course, the energy dependencies 
of soft and hard hadron interaction radii are different, the hard interaction radius 
of hadrons is energy independent, while the soft hadron interaction radius rises 
with energy as In s, due to this difference in the energy dependencies the limit 
Rg(s) — ► at s — > oo takes place. 

After these qualitative remarks we would like to provide a quantitative esti- 
mates for the energy dependence of Rg(s) for DPE processes. We adopt a most 
simple choice and assume that the condensate distributions D^* over impact pa- 
rameter b is controlled by the pion mass, then the convolution 

D h c } <g> D h c 2 ~ exp(-m^). (21) 

In Eq. (1201) we do not write down explicitly dependence on the parton momentum 
fraction x; one can suppose that the x and b dependencies are factorized (at the 
[/-matrix level) and give contribution to the cross-section cr^(s, Q which will be 
cancelled out in the dependence of R 2 s (s). Now we calculate this dependence with 
the models for the function U(s,b) and for the impact parameter dependence of 
the function cth(s, b, £). According to Eq. (T2T|) we should calculate the ratio of 
the two integrals er?/(s, b, ()r](s, b) and o//(s, b, Q over the variables b and The 
function R 2 (s) in the range from the Tevatron energies till the LHC energies and 
beyond this region is depicted on Fig. 3. 

Of course, numerical predictions for the gap survival probability obtained here 
are model dependent, but the qualitative picture of the energy dependence of R 2 (s) 
reflects transition to the new scattering mode at the LHC energies. It is interesting 
that situation is more favorable at the LHC energies since the obtained numer- 
ical values of R 2 (s) at these energies are close to maximal values around 0.2 
and this should lead to higher cross-sections e.g. for Higgs production in DPE 
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Figure 3: Energy dependence of survival probability Rg(s). 



processes compared to the values obtained with eikonal estimations of the gap 
survival probability. Thus, appearance of the reflective scattering would result in 
better detecting of Higgs boson in DPE processes. 

Conclusion 

We would like to note that Eq. © is a more general definition of the gap survival 
probability than Eq. CO). Both are equivalent in the case of absorptive scatter- 
ing, but they are different when the reflective scattering is present. The relevant 
definition in the latter case is Eq. ©. 

The use of Eq. leads to the value of Rg(s) ~ 0.2 for the hard DPE pro- 
cesses at the LHC energies instead of unity which has been obtained in ifTTl on 
the base of the definition Eq. (OQ). The energy dependence of R"i(s) is determined 
by the interplay of soft and hard interaction radii. At s — > oo the ratio R%(s) — > 
due to dominance of the soft hadron interaction radius over the hard interaction 
radius of hadrons. 

The LHC energies appear to lie in the favorable energy region, where soft 
and hard hadron interaction radii are close and this allows one to arrive to rather 
optimistic conclusions on the Higgs finding in the DPE processes there. 
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